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1. Outline 

This is a general introduction into the subject aimed at a general theoretical physics 
audience. Several recent results on the behavior of the cold and hot phases, as well 
as of the critical line have been included, however no attempt was made at being 
exhaustive: further material is to be found in the contributions by A. Nakamura 
and Ph. De Forcrand, as well as in introductory - and topical reviews El. 

This writeup is organized as follows: we give a first introduction and physical 
motivation. Next, we briefly review field theory thermodynamics, and the lattice 
approach. Technical difficulties (the sign problem) are introduced there, while the 
three main methods to circumvent them are discussed in the following Section. 
The following section discusses in some detail the interrelation between the Taylor 
expansion and the imaginary chemical potential approach. Section VI, VII, VIII 
give an overview of the results on the critical line, hot and cold phases respectively. 
A short summary concludes the paper. 

2. Introduction 

The historical developments of the phase diagram of CD is characterized by an 
increasing complication: early views were based on asymptotic freedom , and divided 



February 1, 2008 5:10 WSPC/INSTRUCTION FILE 



mpl 



2 M.P. Lombardo 

sharply the phase diagram into an hadronic phase and a quark gluon plasma phase. 
In the late 90's it has been appreciated that the high density region is much more 
complicated than previously thought ^ 

In the last couple of years, it was the turn of the region above T c to become more 
rich: the survival of bound states above the phase transition brought up the idea of a 
more complicated phase, the strongly interactive quark gluon plasma, characterized 
by a rich particle spectrum of colored particle (free quarks, diquarks, etc,.)™ 

What has the lattice to say in this scenario? Its main task is to provide ab 
initio calculation of the properties of the phase diagram, including, of course, the 
precise location of the critical line, and the equilibrium properties of the different 
phases. One main challenge, of course, is to frame these results into the more general 
context of real time evolution. In this respect, it is worthwhile to remind ourselves 
that equilibrium solution are steady state solution of the dynamical Fokker Planck 
operator; and that lattice calculations will help validating simple models which can 
be studied out of equilibrium: equilibrium studies are a mandatory first step toward 
a full real time understanding of the collisions. 

3. A brief review of field theory thermodynamics 
3.1. Formulation 

Let us remind ourselves how to introduce a chemical potential /i for a conserved 
charge N in the density matrix p in the Grand Canonical formalism, which is the 
one appropriate for a relativistic field theory ^ 



3.2. The Hamiltonian Formalism 

My task in this lecture is to discuss numerical results for the Lagrangian formalism. 
Still, it is very important to mention the Hamiltonian approach: early studies did 
show a great promise, and the recent efforts might well indicate that this is correct 
avenue to treat the phase diagram at finite temperature and density. I refer the 
reader to recent literature on the subject, which contains a full set of references to 
early work El 

3.3. The Lagrangian Formalism 

The path integral representation of the grand partition function Z in the Euclidean 
space gives the temperature as the reciprocal of the imaginary time: 



P = e 



(H-^N)/T 



(1) 

(2) 





(3) 
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with periodic boundary conditions in time for bosons <p{t = 0, x) — <j){t — 1/T,x) 
and antiperiodic for fermions ip(t = 0,a?) = — ip(t = 1/T, x). 

All in all, Z at finite temperature T and density fi is the partition function of a 
statistical system in d+1 dimension, where T is the reciprocal of the imaginary time, 
and n couples to any conserved charge. This representation, which is the starting 
point for a lattice calculation, allows us to deal with thermodynamics and spectrum 
exactly on the same footing. 

The theory is regularized on a space time lattice: a regular four dimensional grid 
with N s points in each space directions, N t points in the imaginary time direction, 
and spacing a. We refer to the very many excellent reviews and textbooks for back- 
ground material on lattice field theory, and we briefly summarize here the specific 
aspects of lattice QCD thermodynamics which will be useful in the following. 

The temperature T on a lattice is the same as in the continuum: T = 1/N t a, 
N t a being the lattice extent in the imaginary time direction (while, ideally, the 
lattice spatial size should be infinite). A lattice realization of a finite density of 
baryons, instead, poses specific problems: the naive discretization of the continuum 
expression fi^jQtp would give an energy e oc ^ diverging in the continuum (a — ► 0) 
limit El. 

The problem could be cured by introducing appropriate counterterms, however 
the analogy between fi and an external field in the 0th (temporal) direction offers 
a nicer solution by considering the appropriate lattice conserved current This 
amounts to the following modification of the fermionic part of the Lagrangian for 
the Othdirection L° F : 

L° F {P) = ^7o^ a VW6 - 4+67oe-^x (4) 

while the remaining part of the Lagrangian is unchanged. This yields the current: 

Jo = -d,L = -d^Lp ^) = 4> xlo e^ x+d + ^ +6 7oe^ a ^ (5) 

This representation of Jo is amenable to a simple interpretation: the time forward 
propagation is enhanced by e^" 1 , while the time backward propagation is discouraged 
by e^"; hence, the link formulation generates a particles-antiparticles asymmetry. 
In addition, note that J Jo = N — N as it should. An alternative way to look at 
the link formulation introduces an explicit dependence on the fugacity e M / T via 
an unitary transformation for the fields El. J n this way L(fi) — L(0), and the /i 
dependence is on the boundaries, via the fugacity e M ^ T : ^(x+Nt) = —e^ aNT ip{x) = 
—e^/ T ip(x). This is analogous to the continuum casJ3 



3.4. Calculational Schemes 

Having set up the formalism, the task is to compute 

Z = [ dUdiPe- s{u ^ 



(6) 
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where from now on the Lagrangian defining the Action will be that of lattice QCD, 
containing gluon fields U and quark fields ip. 

We have two options. We might integrate out gluons first: 

dUdi>dip2(T,n,ip,ilj,U) ~ J dipdipZ(T,fi,<ip,i/)) (7) 

This produce an effective approximate fermion model: the procedure is physically 
appealing, but not systematically improvable, but for one special (lattice) case (see 
below) . Alternatively, we might integrate out fermions exactly, by taking advantage 
of the bilinearity of the fcrmionic part of the Lagrangian L = Lym + Lp = Ly m + 

•4>M(U)ip : 

dUdipd$Z(T,Li t ,i),i>,U)= [ dtf e -(s™([/)-; os (detM)) ^ 



The "effective" model we build this way is exact: the price to pay being that its 
physical interpretation is not as clear as for effective fermion models. Anyway, this 
expression is the starting point for numerical calculations: the fact that in many 
cases they are highly successful tell us that the configuration space is well behaved 
enough that only a minor subset of configurations, although carefully chosen via 
importance sampling, suffice to produce reasonable results. 

3.5. Effective Fermionic Models: analytical approaches 

Let us start by following the first idea, namely integrating out the gluon fields so to 
define an effective fermionic Action. This is a time honored approach, leading, for 
instance, to the instanton model Hamiltonian, hence to the exciting discoveries on 
the QCD phase diagram of the last five years^l . 

On the lattice, one very interesting approach leading to a fermionic model is 
provided by the strong coupling expansion: in the infinite gauge coupling limit the 
Yang Mills term decouples from the Action, and the integral over the gauge fields 
can be carried out exactly. 

The starting point is the QCD lattice Lagrangian: 

3 

S = -l^^^^MlxW^Wxfx+j) - x(x+j)U](x) X (x)] (9) 

X j = l 

-1/2 Vo(x)[x(x)Uo(x)x(x + 0) - X (x + 0)U&) X {x)] 

X 

4 

-l/3]T6/. g 2 J2 [1-reTrU^x)] 

x ^,i/=l 

+ ^2mxx 

X 

The X> x are the staggered fermion fields living on the lattice sites, the C/'s are 
the SU(N C ) gauge connections on the links, the 77's are the lattice Kogut-Susskind 
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counterparts of the Dirac matrices, and the chemical potential is introduced via 
the time link terms e M , e _M as discussed above. This time we have written down 
explicitly the lattice Action to show that the pure gauge term 
Sq = —l/3^:r 6/g 2 X)/i i/=i 4 [l — re TrU^ v (x)\ contains the gauge coupling in the 
denominator, hence it disappears in the infinite coupling limit. Consequently, one 
can perform independent spatial link integrations, leading to 

Z = J Yl dU t d X dxe~ 1/AN ^<^y>^ x]x{xmy]x{v) e- St (10) 

timelinks 

where ^2 <x y> means sum over nearest neighboring links, terms of higher order 
have been dropped, and we recognize a four fermion interactional Further manip- 
ulations yield the mean field effective potential: 

V eff (< tpip >,n) = 2cosh(rN t N c fi)+smh[(N t +l)N c < ^ip >}/sinh(N t < tpip >)(11) 

which we quote for further reference. A standard analysis of V e ff finally gives the 
condensate as a function of temperature and density, and allows the reconstruction 
of the phase diagram. 

More recently this approach has been furthered both in two and three 
color and new developments on cluster algorithms have appeared as well"'. 

In order to describe in detail the rich physics of the finite density phase, one needs 
both to include higher order terms into the strong coupling expansion, as well as to 
go beyond a simple mean field analysis, which assumes an homogeneous background. 
The question is as to whether such improved strong coupling approaches would be 
able to generate a four fermion term with the correct flavor structure as well as order 
of magnitude, thus opening the possibility of a systematically improvable approach 
to finite density QCD, including the study of the superconducting phase. 



3.6. Effective Gluonic Models: Importance Sampling and the 
positivity issue 

Let us write again 

Z(T,a) = J dUe~ iSYM{U) ' lo9idctM)) (12) 

When det M > the functional integral can be evaluated with statistical 
methods, sampling the configurations according to their importance (Sym{U) — 
log(det M)). For this to be possible the would-be-measure (det M) has to be posi- 
tive. 

Let me mention at this point that the factorization 

methodEEIl 

might alleviate 

the problems of complex measures by guiding the simulations along a sensible path 
in the phase space. I will not dwell on this interesting development which is not 
really in the scope of an introductory review, but I wish to call on it the attention 
of the interested reader, as it really seems to offer some promise, and has been 
already tested in random matrix models. 
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In QCD with an even number of flavors, and zero chemical potential, standard 
importance sampling simulations are possible if det M is real, which is true if M* = 
-PMP- 1 where P is any non singular matrix. In the most popular lattice fermion 
formulation this holds: for Wilson fermions P = 75 and for staggered fermions 
P = I (note that this basically expresses a particle-antiparticle symmetry). We will 
consider staggered fermions from now on. 

Consider now the relationship M'(/j,b) — — M(— /is) implying that reality is 
lost when Re/J, 7^ 0: the reality of the determinant is lost, and with it the possibility 
of doing simulations with non zero chemical potential, when we want to create a 
particle antiparticle asymmetry. On the other hand a purely imaginary chemical 
potential does not spoil the reality of the determinant: indeed, even if an imaginary 
chemical potential can be used to extract information at real chemical potential, it 
does not create any real particle-antiparticle asymmetry and it is natural that the 
fermion determinant remains real. 

Note that in QCD with two color the determinant remains positive with nonzero 
real chemical potential: indeed, in that case quarks and antiquarks transform under 
equivalent representation of the color group and are, essentially, the same particle. 
Other important models wit h a real determinant include finite density of isospin^l 
and four fermion models 

All in all, if we want to extract information useful for QCD at nonzero baryon 
density by use of standard MonteCarlo sampling we will have to use information 
from the accessible region: 

Re fi = 0, Im fj, < 
4. Overview of the methods 

To begin with, it is useful to think of the theory in the T, n 2 plane. Let us then 
discuss the phase diagram from the perspective of analyticity and positivity of the 
partition function and of the determinant. One important consideration to keep 
in mind: the Gran Canonical partition function has to be positive. It is only the 
determinant which can change sign, or even be complex, on single configurations. 

Let us consider a mapping from complex [i to complex fi 2 . Because of the sym- 
metry properties of the theory, this mapping can be done without loss of generality. 
Let us note then that Z(p 2 ) is real valued for real fi 2 : this is a situation familiar 
from condensed matter: the partition function is real where the external parameter 
is real, complex otherwise. 

The reality region for the partition function represents states which are physi- 
cally accessible. The reality region for the determinant represents the region which 
is amenable to an importance sampling calculation: Re/i 2 < 0. The methods which 
have been applied so far are 

• [M = Derivatives, Reweighting, Expanded reweighting 

• /1 2 < Imaginary chemical potential 
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Hadron Gas 
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Fig. 1. Sketch of the phase diagram in the /i 2 ,T plane: the solid line is the chiral transition, 
the dashed line is the Roberge Weiss transition. Simulations can be carried out at fi 2 < and 
results continued to the physical domain fi 2 > 0. The derivative and reweighting methods have 
been used so far to extract informations from simulations performed at fi = 0. The imaginary 
chemical potential approaches uses results on the left hand half plane. Different methods could be 
combined to improve the overall performance. 

4.1. Derivatives at fi — 0.0 

This is one early attempt at exploring the physics of nonzero quark density: the 
derivatives can be formally computed at \i = ^! The obvious limitation is that 
we do not really know how far from the /i = axis can we get. Nonetheless, such 
derivatives are interesting per se, and the region where derivatives are clearly dif- 
ferent from zero is the natural candidate for the application of other methods. 

We would also like to quote new results from the Derivatives on the mass spec- 
trum where the first order response of the nucleon (neutron) to baryochemical 
potential was computed in correspondence of several quark masses. 

4.2. Reweighting from \i — 

Back in the 80's Ian Barbour and collaborators proposed to calculate Z(fi) from 
simulations at /i = 0: 



In other words, the chemical potential \i of the target ensemble at that of the 
simulation ensemble - \i = - are different: the properties of the target ensemble 




(13) 
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can be inferred from those of the simulation ensemble, provided that there is a 
sizable overlap between the two ^1 

At T = the Glasgow procedure fails because of a poor overlap (aside, the 
strong coupling calculations were quite useful to asses these problems), and it is 
instructive to study the overlap problem as seen in the Gross Neveu model, where 
there is no sign problem and the results o bta ined with reweighting methods can 
be compared with those of exact simulations 

The distribution of the order parameter (the a particle) helps visualizing the 
problem the order parameter distributions in the two phases do not overlap^. 

The conclusion from these early studies was that reweighting fails in QCD at 
zero temperature because of a poor overlap, and that the reason behind the failure 
is practical rather than conceptual: the situation can be ameliorated if a better 
starting point were used. 



4.3. Fodor and Katz's multiparameter reweighting 

The prescription for ameliorating the overlap is due to Fodor and Katz whose 
Multiparameter reweighting use fluctuations around T c at \x — to explore the 
critical region. Making reference to Fig. 2, and oversimplifying: instead of trying to 
reweight the distribution at zero temperature in the broken phase, which is obviously 
hopeless, one might hope that a distribution generated at zero density, and close to 
the critical temperature, bears more resemblance with the target distribution along 
the critical line, and is thus amenab le to a successful reweighting. 

The strategy was applied to QCDI232II. The improvement obtained is impressive 
and produced the first quantitative results for the critical line at nonzero chemical 
potential in QCD: we will come back to this in the section on results. A multistep 
reweighting proposed by Crompton might well produce a further improvement. 



4.4. Taylor Expanded Reweighting 

The Bielefeld-Swansea collaboration suggested a Taylor expansion of the reweighting 
factor as a power series in A = /i/T, and similarly for any operatoJ^EHl _ 

This strategy is computationally very convenient as it greatly simplifies the 
calculation of the determinant. Expectation values are then given by 



, , _ ((Oo + OiX + O2X 2 + ■ ■ ■) exp^A + K 2 \ 2 + ...- AS g )) x=0iPo 
{ HM <cxp(7e 1 A + ft 2 A 2 + ...-AS s ))A=o,/ 3 o 



Results have been obtained both for the critical line and thermodynamics. 



4.5. Imaginary Baryon Chemical Potential 

This method uses information from all of the negative /s 2 half plane (Fig. 1) to 
explore the positive, physical relevant region. An imaginary chemical potential v in 
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a sense bridges Canonical and Grand Canonical ensemble^! 

ZcW = 7T / dvZ GC (i V )e-^ vN (15) 
27r Jo 

The main physical idea behind any practical application is that at /i = fluctuations 
allow the exploration of Nb ^ hence tell us about \x ^ 0. Mutatis mutandis, this 
is the same condition for the reweighting methods to be effective: the physics of the 
simulation ensemble has to overlap with that of the target ensemble. 

A practical way to use the results obtained at negative u ? re lies on their an- 
alytical continuation in the real plane. For this to be effective^' 2(fi, T) must be 
analytical, nontrivial, and fulfilling this rule of thumb: 



X (T, /i) = dpfji, T)/dn = d 2 logZ{^ T)/V > (16) 

las I 
:ed i 

the two color model 



This approach has been tested in the strong couplin g lim it of QCD, in the 
dimensionally reduc ed m o del of high temperature QCD I2HI and, more recently, in 



5. Taylor expansion and analytic continuation 

We discuss here in more detail the interrelation between the Taylor expansion and 
the imaginary chemical potential approach. 

Res ults o btained at imaginary fiB can be analytically continued to real 
/zJ2E32E2. i n principle, rigorous arguments guarantee that the analytic continu- 
ation of a function can be done within the entire analytic domain. In practice, the 
exact analytic form is not known, and a systematic procedure relying on the Taylor 
expansion is only valid within the circle of convergence of the series itself. Here, 
we discuss how to implement the analytic continuation of the critical line and of 
thermodynamics observables beyond the circle of convergence of the Taylor series 
in a controlled way. 

Let us remind ourselves that an analytic function is locally representable as a 
Taylor series. The convergence disks can be chosen is such a way that they overlap 
two by two, and cover the analytic domain. Thus, one way to build the analytic con- 
tinuation is by connecting all of these convergence disks. The arcs of the convergence 
circles which are within the region where f is analytic have a pure geometric mean- 
ing, and by no means are an obstacle to the analytic continuation. Assume now that 
the circle of convergence about z = (0,0) has radius unit, i.e. is tangent to the lines 
which limit the analytic domain; take now a z value, say z\ = (0,a), 1/2 < a < 1 
inside the convergence disk as the origin of a new scries expansion, which is explic- 
itly defined by the rearrangement (z — zq) u — (z — Z\ + z\ — zo) n As the radius 
of convergence of the new series will be again one, this procedure will extend the 
domain of definition of our original function (the two series define restrictions of 
the same function to the intersection between the two disks), and by 'sliding' the 
convergence disk we can cover all the analytic strip. 
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We have skecthed above the standard theoretical argument to demonstrate the 
feasibility of analytic continuation beyond the radius of convergence, and we will 
show that the Pade' series is one practical way to accomplish it. 

To complete this discussion, let me mention that the radius of convergence of 
a Taylor expansion about the origin might well be larger than the distance of the 
origin itself from the nearest singularity. While complex analysis textbooks offer a 
full discussion of this point, I would like to single out here three cases which might 
be encountered in usual critical behavior: fl(z) = ^4l(z)(l — z/z c )~ x ; f2(z) = 
A2(z)9(z — z c ); f3(z) — Ai(z)9(z — z c )(l — z/z*)~ x , where An(z) is an analytic 
function. 

Case 1 corresponds to an usual critical behavior (second order or larger). Case 2 
represents a strong first order phase transition. Case 3 is intermediate between the 
two, a weak first order transition at z Cl and a spinodal point at z* . a Correspondingly, 
we have different radius of convergence of the Taylor series: ri = \z c \, T2 = oo, r% = 
z* : in conclusion the radius of convergence of the Taylor expansion for the critical 
line and thermodynamics observables might be infinite as well a finite, depending 
on the nature of the Roberge Weiss transition. Conversely, if the nature of the phase 
transition is known, one can infer from it the radius of convergence of the Taylor 
series, as done by Gavai and GuptaE^, which in turn locates the critical point. 

6. The Critical Line : Taylor vs Pade' 

A full discussion of the current status of the critical line we refer to recent reviews^. 
Here, we would like to emphasize more general aspects on the possibility of contin- 
uing the results beyond the radius of convergence of the Taylor expansio, by taking 
adantage of the discussions preented in the Section above. 

The radius of convergence of the Taylor representation of the critical line might 
well be limited by the Roberge Weiss singularities (see again Fig.l). However, as 
explained before, the Pade' appoximation is not. 

WeEnih ave performed present the Pade' analysis where we have used data for 
four 021 and two flavor The results seem stable beyond hb — 500MeV((J,B /T ~ 
1), with the Pade' analysis in good agreement with Taylor expansion for smaller 
/i values. At larger [i the Taylor expansion seems less stable, while the Pade' still 
converges, giving a slope of the critical line larger than the naive continuation of 
the second order Taylor approximations. The same behavior is sug gested by recent 
results within the canonical approacr!^ and the DOS method 

We underscore that the possibility of analytically continue the results beyond 
the radius of convergence of the Taylor series by no means imply that one can 
blindly extrapolate a lower order approximation! Even when it is possible to achieve 

a The analytic continuation is insensitive to a discontinuous phase transition since it lives on 
the metastable branch : it follows the secondary minimum and determines the spinodal point 
< "ipip >= A([i — /i*)' 3 . The discontinuity can be related to fi — fi* . Both shrinks to zero at the 
endpoint of a first order transition. 
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convergence - via Pade' approximants, or within the convergence radius of the 
Taylor series - one has always to cross check different orders of approximation to 
make sure that convergence has indeed been achieved. For instance, it would then be 
interesting to repeat the comparisons between the second order results shown inEH 
by extending the Taylor series to fourth order and/or by use of Pade' approximants. 



6.1. New Results for Wilson Fermions 

It is very important to double check the results coming from different meth- 
ods and different approaches. The results decsribed so fare have been ob- 
tained with staggered fermions. First calculations with Wilson fermions at 
fintie density, and imaginary chemical potential, have appeared recently- 
Chen:2004tb,Luo:2004se,Luo:2004mc. The results are so far in nice agreement with 
those obtained with the staggered formulation. 



7. The Hot Phase and the approach to a Free Gas 

The behaviour of the number density at high temperature approaches the lattice 
Stephan-Boltzmann prediction, with some residual deviation. The deviation from a 
free field behavior can be parametrised as as^IESI 

AP(T )jU ) = /(T, M )P/ ree (T,M) (17) 

where Pf ree {T, /i) is the lattice free result for the pressure. For instance, in the 
discussion of Ref. GUI 

f(T,u) =2(l-2a,/7r) (18) 

and the crucial point was that a s is fi dependent. 

We can search for such a non trivial prefactor f(T, /i) by taking the ratio between 
the numerical data and the lattice free field result n^ ree (^i) at imaginary chemical 
potential: 

R{T,a I ) = ^f\ (19) 

A non-trivial (i.e. not a constant) i?(T, uj) would indicate a non-trivial f(T,/ij. In 
Ref. we calculated R(T,/i]) versus Hi/T: the results for T > 1.5T C seem consistent 
with a free lattice gas, with an fixed effective number of flavors N'f* '(T) / '4 = R(T): 
Nf f = 0.92 x 4 for T = 3.5T C , and N e / f = 0.89 x 4 for T = 1.5T C . 



7.1. Beyond fi/T ~ 1 in the Hot Phase 

The Pade' approximants to the results for the chiral condensate in the hot phase 
have been calculated as well^l using four flavor data^l Again we found that the 
Pade' analysis seems capable to produce stable results. We should also note that 
the Taylor expansion seems stable as well, which might indicate a large (infinite?) 
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radius of convergence in this range of temperature. Indeed, as noted in 1^ the radius 
of convergence should tend to infinite in the infinite temperature limit , an d indeed it 
has been estimated to be large by the Bielefeld-Swansea collaboratiorP-4 A detailed 
investigation at imaginary /i of the region closer to the critical temperature is in 
progress t±al . 

We also note a possible interplay of thermodynamics and critical behaviour for 
Tc < T < Te — 1.1T C : the critical line at negative fi 2 would imply , at least for 
second order and weak first order transitions, logP(fi, T) oc (/i 2 — n 2 ) 11 which is 
incompatible with a free field behaviour. 

8. The Hadronic Phase the Hadron Resonance Gas Model 

The Hadron Resonance Gas model might provide as description of QCD thermo- 
dynamics in the confined, hadronic phase of QCD: the grand canonical partition 
function of the Hadron Resonance Gas model"-^ has a simple hyperbolic cosine 
behaviour. 

This behaviour could be assessed via a computation of the Taylor coefficients 1*^1 
As a perhaps simpler alternative, the hadron gas model can be framed in the discus- 
sion of the phase diagram in the temperature-imaginary chemical potential plane 
which suggests to use Fourier analysis in this region, as obscrvables arc periodic and 
continuous thcr 

For observables which are even (O e ) or odd (O ) under [i — > — /i the analytic 
continuation to real chemical potential of the Fourier series read O e [o](fxj, N t ) = 
yZr, a F t " > cosh[sinh|(niVt./V e /x/). In a Fourier analysis of the chiral condensate 1^ 
and of the number density^ - even and odd observables, respectively - we lim- 
ited ourselves to n — 0, 1,2 and we assessed the validity of the fits via both the 
value of the x 2 /d.o.f. and the stability of dp and given by one and two 
cosine [sine] fits: when HRG holds true, one term in the Fourier series should 
suffice. (sinh(x) — > sin(x)) n(n) = — Ksin(N c N t ^). 

We found that one cosine [sine] fit describes reasonably well the data up to 
T ~ 0.985T C ; further terms in the expansion did not modify much the value of 
the first coefficients and does not particularly improve the x 2 /d.o.f.. This means 
that our data are well approximated by the hadron resonance gas prediction AP cx 
(cosh(^ B /T) - 1) in the broken phase up to T ~ 0.985T C . 

In the same region, we can also computer the mis ma tch with respect to the 
HRG in an 'effective mass analysis' style: this analysis EH confirms that the HRG 
is consistent with our data within errors up to T — .985T C 

Mismatch = n(n) / 1 sin{N c N t ^) — k (20) 
8.1. The critical line from the Hadron Gas 

An alternative way to analytically continue the results relies on phenomenological 
modeling. The Hadron Resonance Gas model might provide a description of QCD 
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thermodynamics in the confined, hadronic phase of QCD "" * ' I l> and can be used 
to determine the critical line as well. 

The critical temperature as a funct ion of \ib is determined by lines of constant 
energy density: e ~ 0.5 — 1.0 GeV/fm^S3 A continuation of the critical line using 
the HRG ansatz plus a fixed energy (or any other quantity determined at fj, = 
) criterion suggests the implicit form for the critical line T = f(T)cosh(pB/T) 
with lim^/7^0 f (T)cosh(nB /T) = 1 — fc/i 2 . We have naively approximated f(T) = 
1 — kfi 2 , and used the resulting form to fit the data in the fi/B < 1 range. Accrdoing 
to the above discussion, this again can be continued beyond this limi t, a nd also in 
this case we get a critical line whose slope increases with increasing /JsJa. 



8.2. How to calculate the critical values 

These calculations described above require the critical density as an input. Let us 
then consider: 

n(iu) — aisin(ifiN c NT) + a2sin(i2fiN c Nt) (21) 

Analytic continuation up to fi = ji c (T)) gives: 

n{n) = aisinh(^N c N t ) + a 2 sinh(i2fiN c N T ) (22) 

This gives the critical density at T — .985T C and mass = .05 in lattice units: 
n c (fj, c )/T 3 ~ o,5 IMI once the value of fi c has been taken into account. In addition, 
the mass dependence of the critical density has been be extimated from the 
Maxwell Relations: 

d < > /du = dn{u)/dm (23) 
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9. Parting Comments 

We have described three different, independent methods which afford a quantita- 
tive study of the phases, and phase transitions in QCD. All of the methods exploit 
physical fluctuations either at zero and purely imaginary chemical potential to ex- 
plore real baryon chemical potential. The physical idea in a sense is similar, but the 
systematics is very different. Cross check are then most useful and informative, and 
in many cases, which we have reviewd, have been performed satisfactorily. 

I should reiterate a caveat: both the critical line, and thermodynamics observ- 
ables in the hadronic phase are very sensitive to the quark masses. I preferred to 
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concentrate on methods and general idea, rather than on rapidly evolving results. I 
apologize for possible incomplete presentation, and refer once more to the reviews I 
cited. The reader who chances into this note at some later stage, might find useful 
to read the plenary reviews on Thermodynamics at the Latticexx Conferences and 
that on Lattice at the QuarkMatterxx Conferences. 

Finally, again and again, the methods described here are dodges for the sign 
problem reviewed above. A more complete solution might well be afforded by the 
Hamilt onian a pproach reviewed™ the density of sta tes formalism EH, the canonical 
method 0231111 or the strong coupling expansion 
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